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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO . 1187 


EQUATIONS FOR ADIABATIC BUT ROTATIONAL STEADY 
GAS FLOWS WITHOUT FRICTION* 

By Manfred Sch&fer 


ASSUMPTIONS 

1. The flowing gases are assumed to have uniform energy distri- 
bution. ( "isoenergetic gas flows," that is, the energy equation 

is valid with the same constants for the entire flow.) This is 
correct, for example, for gas flows issuing from a region of con- 
stant pressure, density, temperature, and velocity. This property 
is not destroyed by compression Bhocks because of the universal 
validity of the energy law. 

2. The gas behaves adiabatically, not during the compression 
shock itself but both before and after the shock. However, the 
adiabatic equation (p/p K = C) is not valid for the entire gas 
flow with the same constant C but rather with an appropriate 
individual constant for each portion of the gas. For steady flows, 
this means that the constant C of the adiabatic equation is a 
function of the stream function. Consequently, a gas that has 
been flowing "isentropically", that is, with the sane constant C 
of the adiabatic equation throughout (for example, m origination 
from a region of constant density, temperature, and velocity) no 
longer remains isentropic after a compression shock if the com- 
pression shock is not extremely simple (wedge shaped in a two- 
dimensional flow or cone shaped in a rotationally symmetrical flow) . 

The solution of nonisentropic flows is therefore an urgent 
necessity. 

For the computations, the physical system of measurements 
(mass, length, time, degrees Kelvin (absolute temperature), and 
calorie) will be used. This system was chosen in order to avoid 
the frequent occurrence of g. In the final equations, the gas 
condition no longer explicitly appears but instead only properties 
of the velocity field based on the "limit velocity" V m (corre- 
sponding to outflow into a vacuum) appear. 


*"Gleichungen fur Adiabatische, aber Wirbelbehaftete Stationare 
Gas stri'mungen ohne Reibung." Lehrstuhl fur Technische Mechanik an 
der Technischen Hochschule Dresden, Archiv. Nr. 44/l. 
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SYMBOLS 

p pressure 
p dens ity 

T absolute temperature 

v velocity vector [MCA comment: Changed from the German script 

v for simplicity.] with velocity components u, v, and w 
in a cartesian x, y, z coordinate system 

q magnitude of velocity 

i | i stream function 

limit velocity corresponding to outflow into a vacuum 
a local sonic velocity 


Motion equations: 


BASIC EQUATIONS 



5u X| , Bp 

Bz J 3x 



( 1 ) 


p v • V v = - grad p 
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Continuity equation: 

<*> 

or 

div (p v) =0 

Along each stream line \|/ = constant, the following adiabatic 
equation applies: 

Adiabatic equation: 

- fW) (3) 

P 

The flow line has direction cosines that are related by the propor- 
tion u : v :wj hence 

^u + ^-v + ^w = grad • v = 0 (4) 

dx dydz 

However, the stream function is not completely determined by this 
property. The numbering of the stream lines will be on the basis 
of the mass flowing through, in the cases of two-dimensional and 
of rotationally symmetrical flows. 

The following gas equation is understood to apply: 

Gas equation: 


£ * RT (5) 

P 

The quantity or, correspondingly, f (ty) is related in a 

readily apparent manner to the potential temperature used by 
meteorologists as well as to the entropy, namely: 

The entropy per unit mass is 


S 


1 B 
J K - 1 


In 
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J, the mechanical heat equivalent, equals 


4.ie x io 7 6 ram . St -22S± 

cal 

The "potential temperature" 9 is that temperature which a 
gas would reach if it were adiabatically brought to a suitably 
selected normal pressure p n . 

In accordance with the adiabatic equation 

JL In _ !il i-« 

P K ' ' P/ P ” 

In accordance with the gas equation 


Jj,. 8***1-* 

Consequently 

„K Pn^ 1 P 
" E k p k 

Instead of the potential temperature, F . Frankl introduces a 
"zero density" p 0 . (See discussion of the literature at end of 
report . ) 

The energy equation per unit mass flow is, as proved in thermo- 
dynamics, 


JE.^ + £ + i (u 2 + v 2 + w 2 ) = constant 

where, according to the introductory remarks, the constant is supposed 
to be valid for the entire gas. The internal energy in thermal units 

is E-i ; je- + E is the heat content i per unit mass, expressed in 
1 P 

mechanical units. 
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In the ideal gases 


T? — 

B i “ 


c v T 


c v 


1 _R_ 
J K-l 


i = JEj + E 
1 P 


= JRT + E 

K-l p 



1 

K-l 



- E Ji_ 

p K-l 


Hence the energy equation becomes 


^ i(u 2 + v 2 + w 2 ^ = constant 


In terms of the maximum possible velocity V m , which is the 

V 2 

"limit velocity," -2— may be written for the constant; for with 

Lt 

outflow into a vacuum (p/p = 0), according to the previous 
equation 



+ v 



Y 2 
v m 


2 


Thus the final form of the energy equation is obtained if one more 
quantity, the magnitude of the velocity q is introduced. 

2 V 2 

£ J5_ + £_ = _HL (6) 

p K-l 2 2 

n.b. For adiabatic flows, the following is true along a flow line; 



= f«r)K 



K-2 


dp = 


K 

K-l 


ffo)p K ‘ 


JL. E 

K-l p 
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and this equation agrees with the first term of equation (6). 
Because J' ^ occurs in the integration of the motion equation (1), 

v • V v it is evident that the energy equation (6) may 

he derived hy this method also, at least for adiabatic flows. 
Compared with this, however, the thermodynamic derivation used 
previously has the advantage of being valid for the entire flow 
independently of compression shocks. 

By the use of equation (6), the following equation is obtained 
for local sound velocity a: 

a 2 = g - K fMp*- 1 * K ^ (v m 2 - q 2 ) (7) 


DIFFERENTIAL EQUATIONS FOR GAS-VELOCITY FIELD ALONE 

From equations (l), (2), (3), and (6) it is now intended to 
derive equations that contain only velocities; that is, p and p 
are to be eliminated. 

/ 

In order to eliminate p from equation (2), div(pv) = 0, 
the following equation is derived from equations (3) and (6): 


therefore 


and 


f(\|/)p 


K-l 


K 

K-l 



P 


I JSii 
L 2 K 


v m 2 - <i 2 

fWO 


i 

k-i 


div(p v ) = p div v + v • grad p = 0 


(2a) 


This equation may be replaced by the equivalent equation 
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Now with the new expression for 
tion in this section . 3 

grad In p = grad In [v m 2 
Then it follows that 


p [NACA comment: The second eqaa 

1 1 
- q 2 ] 1 - grad In fftO*” 1 


v • grad In p = v • grad In [v m 2 



because 


x 

crad In fNO*" 1 = -ir df>W — lu 2Sk + v ^ + w 
grad m XW) = ^ + V 9 y + W g z J 


and according to equation (4) becomes equal to 0. Therefore 
div v + v • grad In p = div v + v • grad In [v m 2 - q 2 ] 


1 

K-l 


_ v • grad [V m 2 - q 2 ] 

= div v + = 0 


p p K -1 

C T „, 3 


For this can finally be written 


1 1 
P p K-l _ _ _ p p K-l 

[V m 2 - q^] div v + v • grad [V “ - q d ] =0 


or: 


4iV (v [V m 2 - q 2 ]""^ = 0 

This equation is already found in the wort of L. Crocco, 1937. 
(See discussion of literature at end of report.) 


( 8 ) 
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From the motion equation 


p v * V v = -grad p 


p and p are also eliminated. 
From previous reasoning 


k- 1 ~^m' 
. 2k 


2 - q 2f - 1 


According to equation (3) 


P = f«,) K - 1 |f^ (V„ 2 - q 2 ) K_1 


Therefore — £ is expressed as: 


grad p = grad 


jfOl/f (V m 2 - q 2 ) 


= - th fW/) K_1 ( V • * 2) _ K_1 grad * 


' K_1 W 2 ~ <1 2 )_ K_1 IF grad q* 

K K 

= - jrj ft*)’ 1531 [fr < v n 2 - « 2 >J k ' 1 8 rad * 


ffo) K_1 (V „ 2 - qZ^^grad £ 
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= (v m 2 - <1 2 ) erad^ + grad f 

Consequently, after taking into account also 

2 

v - V v = grad + rotation v X v 
the following is obtained as the motion equation: 

rotation v X v = ^ (V m 2 - q 2 ) — -- grad \|/ (9) 

This equation appears to be new. It contains the expression 
of the Bjerknes vortex theorem for the flows in question. 

Because 

grad ^ - grad In f(40 

L. Crocco's (1937) equation number (5) may also be immediately 
derived from this equation 


rotat ion 




however, use of this equation is less informative than equation (9) 
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INTRODUCTION OF DIMENSIONLESS MAGNITUDES 

If the procedure of other authors Is followed (Castagna: Atti 

della R. Acc. delle Science di Torino, vol. 70) in basing the veloc- 
ities on which appears to be the most useful procedure, the 

following system of equations is finally obtained 


rotation v x v = ~ (1 - q 2 ) grad In f(\|/) 

C ft 


divV v [1 - q 2 ] J = 0 


[NACA comment: Compare equation (8).] and equation (4) as the 

third 


v • grad \J/ = 0 


(III) 


These equations furnish five scalar equations for the four 
unknown quantities u, v, w, and ''V. This apparent extra relation 
is caused by the fact that use has already been made of an integral 
of the gas-dynamic equations, namely the energy integral, equation (6). 
For this reason, the apparent extra relation introduces no contra- 
dictions and proves to be no impediment in practice. 

For brevity, equation (7) is also used in dimensionless form 


a 2 . ^ (1 - q 2 ) 


APPLICATION OF EQUATIONS (i), (II), AND (ill) TO 
TWO-DIMENSIONAL AND ROTATIONALLY 


SYMMETRICAL FLOWS 


1. Two-Dimensional Flows 


Equations (II) and (III) are satisfied by the expressions 

fu (1 - q 2 )*" 1 = i!/ 
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The vector equation (l) yields the following two equations 
Ydl d In f M 

\^Bx 5y/ — 2 k ' ^ 3x 

and 

fai . aA u „ i (i . .2) a ln | i 

V^Bx By/- 2K ^ ^ ' d¥ Sy 

By expression of the underlined velocities in terms of ty, 
the following single equation is obtained 

K 

fa . s«V i- (1 - q 2 ) K ' 1 d ( 11 ) 

V^Bx By ) 2K H 

and the apparent extra relation no longer exists. 

This equation agrees with Crocco's equation, equation (10), 
which that author could reach only by way of integration processes, 
for the entropy in mechanical units is 


S = 


JL 

K-l 


In 


pff = 


B 

K-l 


In f(\|/) 


Equation (11) now must be transformed, 
equation (10) 


In accordance with 


( 1 - a 2 ) 

Bx v xx * ‘ 


1 

K-l 


- — (1 
K-l V 


q 2 ) 


1 

K-l 


-1 


5q2 

Bx 


By use of equation (7), from the previous equation the following 
is obtained: 


- (1 - 1' 


2jK-1 Qv 


_ = \1/ 
Bx xx 


+ i*. a Ij/ _ V (1 - q 2 ) 


1 

K-l 


2a 2 Bx 


xx 


2a 2 


Bq 2 

Bx 
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Correspondingly, it is true that 


Su 

3y 





h. 

K-l 


(1 




and 


1 

(! - q 2 ) K_1 |^ - + 


1 

h. 3ai , * * u (1 - c 2 /- 1 3j£ 

2 a 2 3 y yy 2 a 2 3y 


In order to transform and a^L-, the starting point is 

taken as either 

„ JL 

Q 2 = u 2 + v 2 = (1 - q 2 ) ^ 1 (\|; x 2 + \|/ y 2 ) 

or 

2 

q 2 (l - q 2 )^ 1 = \j; x 2 + \j/ y 2 

When differentiated with respect to • x, this yields 




2 

K-l 


-1 




_2 

K-l 


4*> <1--^- 




- 2(1 
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Consequently 




2 ". *„2 


- 2 {- 7 *xx + U * X y} 


Corre spondLingly 


d . q 2 ) "' 1 ( 1 . g.ad , 2 {. + u ^j. 


The four underlined equations may be summarized as follows: 
1 




v \cx - 2uv^|/ iy + u 2 \l/ yy 


a 2 -q 2 


(a 2 - u 2 )\l/ xx - 2uv\j/ xv + (a 2 - v 2 )\ji 


xy 


r yy 


a 2 - q 2 


Therefore equation (11) becomes 


5v 3u 
3x Sy 


(a 2 - u 2 )\l/ xx - 2uv\l/ xy + (a 2 - v 2 >|/ 


yy 


(a 2 - q 2 ) (1 - q 2 ) 


i (1 

ZK ' 


„2' 1 K-1 d In f (\J/) 
q > d^ 


This further yields 






= . i (i-£Vi_o 2 r 1 zj njm 


2K 


d\J / 
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and finally 



Kfl 




( 12 ) 


This is the differential equation for the two-dimensional velocity 
field, with which also belong the defining equations, (10) and (7), 
for u, v, and a. All velocities are based on the limit velocity 


In the familiar special case wherein f(\|/) is constant through- 
out the whole flow 


( X ' a 2 )^ ‘ 2 a 2 ^ V + ( X ' " ° (13) 

Because according to equation (i) the flow is then free from 
rotation, and hence with the velocity potential cp 

u = 92 v = ( 14 ) 

dx' 3y 

the following equation is obtained from equation (ll): 
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Equation (15) for cp is analogous in form to equation (13) 
for however, it muBt he remembered in this connection that u 

and v are expressed altogether differently in terms of cp than 
in terms of \J/ (compare equations (10) and (14)). 

Holenbroek and Tschapligin (Hodographenmethode der Gasdynamik) 
have done work involving the comparison of the two differential 
equations (13) and (15). 


2. Rotationally Symmetrical Flows 


The cylindrical coordinates r, ( 0 , and z will now be used. 
The velocity components u || z (axial) and v Jj r (radial) are 
independent of 0) . 

Equation (II) becomes: 


JL 

3z 


u(l 


t K-l 


, '”J + ?sL rT(1 -«* ) 


K-l 


= 0 


where 


q 2 = u 2 + v 2 


This equation is satisfied by: 


1 

ru(l - a 2 )"" 1 = |4 
Sr 


< 


rv(l - q 2 ) 


1 

K-l 


(16) 


dvj ]f_ 

3z 


where \|/(r, z) is the stream function, (compare Stokos's struam 
function for incompressible flows) . 


At the same time, equation (ill) is also satisfied by this 
expression, therefore only the vector equation (i) remains to be 
dealt with. Equation (i) yields: 
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/av _ Bu V _ _1 _ 2 ] * In fflQ 

\^9z dr J — 2tt ^ J Aty 9r 

_ 1 3v _ 9u^ _ _1_ f. _ 2 \ A ln f(^) 3^ 

y9z 9ry- ~ 2tt ^ J d\f 9z 

By expressing u and v in terms of \J/, the following single 
equation is obtained: 

K 

(Bv _ 9_u j _ r_ ( 1 _ 2 V " 1 d In f(M0 fl7 > 

\9z dr 2K \ 4 ) d\)/ ; 


This equation is equivalent to L. Crocco's equation, equation (13), 
which, however, he could only reach by integration. 


The necessity of suitably transforming this last differential 
equation, equation (17), remains. 

From equation (16) is obtained 


r If- - - (1 - K -I - (1 ■ qZ) K - 1 




9 


"i d£ 

d z 


and making use of equation (7) 


- (1 - r U = + 


_ l!f _ vr(l - q 2 ) K " 1 3q 2 


3z” v zz + 2 ^zz” 

2a^ 9z 


2a £ 


By analogy 


9u 

9r 


= M, rr (l- 9 2f (i- q 2)- ^ - 1 (x-a 2 )' K - 1 
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and 


* i - v 4*& - v ■ 4 f 


2 2 

In order to transform and , the starting point is 

taken as 9z 3r 


i 2 = u 2 + v 2 = -l 1 — -ail 
1 ~2 


2 

K-l 


(v)/ z 2 +^ r 2 ) 


or 


q 2 r 2 (1 - q 2 )*" 1 = \|/ z 2 + \J/ r 2 
By differentiation with respect to z 

2 


{(■ - -r 


,2 - 2a!^. ( l - q‘ 


K-i " n *„2 


J 3z 




q'^J 3 z 


- \ z * *r 


X 

= z(l - q^f' 1 r|- V\|/„ + 
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From these equations it follows that 


(l - r{l - = z{- ,* 2Z ♦ u* r2 } 


and finally 


(1 - q 2 ) K 1 Qq 2 ~ ^zz + U ^rz 


2a^ 3z 


a 2 - q 2 


By differentiation with respect to r 


i -asf ' 1 r 2 - 2 an: 


2 r 2 ( X ~1 3q 2 rtf rf'f ~1 

_ V 1_ V J 


K-l 




\J/ \i / + \|/ \J/ 

z rz r rr 


= 2 ( 1 -^) K - 1 r{- 


wl/ + u4/ 
Y rz Y rr 


} 


1 - ^ < x - - 


\ 

2l - vtj/ + u\|/ } - 2 

rz rr/ 


(i - ^ 


(1 - 1 2 ) K X r aq g 
2a 2 3r 


(- vii/ rz + mj/ rr ) - (1 - q 2 ) 


K-l 


2 2 
a^ - q^ 
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From the four underlined equations it follows that 


1 



By use of this expression, differential equation (17) for ^ 
now becomes 

. 1 . 

(a 2 - u 2 )^ zz - 2uv^ rz + (a 2 - v 2 )^ rr - a 2 


or finally 


- - £ (a? - q2) (1 - qa)"' 1 "- 1 Uatiai 

ZK * dty 


(l _ + (x . £\, . *r 

V aV “ a 2 V a 2 / ™ r 
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This equation is the differential equation for the rotationally 
symmetrical velocity field, with which belong also the defining 
equations (7) and (16) for u, v, and a. All velocities are based 
on the limit velocity V m . 

In the familiar special case in which fOJ/) is constant through- 
out the whole flow, 



Because, according to equation (i), the flow is then free of 
rotation and with the velocity potential defined as follows: 


u 


v = 
3z 


( 20 ) 


the following is obtained by insertion in equation (il) (the continuity 
equation) : 


3u 3v 
3z 3r 


+ 






- 2u2 

9z 


2uv . 
3z 


2uv 


du 

3r 


2v 2 


3v 

3r 


= 0 




zz 


cp rr + 



~ u2c P ZZ + 2uvC Prz + y \v = 0 


or 



( 21 ) 


Differential equation (21) for cp is the same, except for the 
3ign of the last term, as that for \|/, equation (19). However, u 
and v are expressed altogether differently in terms of \J/ than in 
terms of cp (compare equations (16) and (20)). 
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THE LITERATURE 

1. Frankl, F. and Aleks jeva, R.: Two Boundary -Value Problems 

Involving the Theory of Hyperbolic Differential Equations of the 
Second Order with Applications to Gas Flows of Supersonic Velocity 
(Russian). Matematiceski Sbornik, 1934, p. 483. There exists a 
German script translation of this work made at that time at my 
request by Mr. Bl. Dolaptschiew of Sofia. 

In this report a differential equation occurs for \|/ in the 
two-dimensional case of the general flow (f(^) £ constant) that 
is copied from a paper by F. Frankl in the collection "Soviet Union 
Rockets" which is unavailable . The differential equation set forth 
by Mr. F. Frankl is aB follows: 


A u 2 k-1 V 2 


4 uv S 2 \[/ + / 

\.± 



\ a* 2 k+1 a* 2 

/ 8x 2 

k+1 a* 2 9x3y ' 

V a»2 

k+1 a* 2 / 

/ay 2 


2 


+ *±i a *2 

2k 






k-1 _w£\ k_1 p ^o 
k+1 a* 2 / ° d\|/ 


In the report from which the quotation is taken, the only statement 
made regarding definitions is that the "zero density" p Q is the 
density obtained by the adiabatic reduction of the velocity to zero. 
Aside from the unprofitable introduction of p Q , F. Frankl' s equa- 
tion may be immediately obtained from equation (18) if the sonic 
velocity a is expressed in terms of the critical velocity or 
critical sonic velocity a* (u 2 + v 2 = a* 2 ) in accordance with 


a 2 = ^±i a* 2 - (u 2 + v 2 ) 

2 2 

The meaning of the other symbols insofar as they differ from these 
used here is obviously as follows : 

w 2 corresponds to q 2 
k corresponds to n . 
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2. Crocco, L.: Una nuova funzione cLi correnta per lo studio 

del moto rotazionale del gas, Eend. della B. Accad. dei Bincei, 

Vol* XXIII, serie 6 a , 1° sem. fasc. 2, 1936. 

Translated as "Eine neue Stromfunktion fur die Erforschung der 
Bewegung der Gase mit Botation", ZaMM 17, p. 1, 1937. 

Crocco's final equations (10 11 ) and 0-3’) agree with equations 
(12) and (18). Mention has previously been made of the differences 
in derivation. 


NOTES 

1. With regard to page 3 it should be noted that up to the 
present time a stream function \J/ has been defined only in the two- 
dimensional and rotationally symmetrical flows, which cases are, to 
be sure, of principal interest. In the general three-dimensional 
case, ^ = constant is to be considered as the equation applying 

to those flow surfaces (built up out of flow lines) on which the 
function f(\{/) is constant. 

The Crocco stream function is used chiefly as the stream function 
for two-dimensional and rotationally symmetrical flows. The numbering 
of the flow lineB is, with the exception of the isentropic case 
(chapters III and IV), not done by Crocco on the basis of the mass 
flowing through, but corresponds instead to his form of the continuity 
equation (8). Seasons of expediency led, even in the case of the 
nonisentropic rotational flows (chapter VI), to reversion from the 
Crocco stream function to the customary stream function. 

2. Eegarding page 6 it is to be noted that the relation between 

the new vortex theorem as expressed in equations (9) and (i) and 
the Bjerknes vortex theorem has recently received proper clarifica- 
tion for the first time in a paper in the Wieselsberger-Gedachtnisheft, 
ZWB (W. Tollmien: Ein Wirbelsatz f&r stationare isoenergetische 

Gasstri'mungen. ) In the same paper the apparent excessive determina- 
tion of the equation system (i), (II), and (ill) is decisively 
explained by the fact that the vortex theorem (l), contrary to first 
appearances, yields only two scalar equations. 

3. The relation between pressure, density, and velocity for 
isentropic flows that will be used frequently in the following 
lectures is appended here. When using dimensioned quantities, in 
accordance with the energy equation (6) 
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K-l p 2 

and in accordance with the adiabatic equation for iaentropic flows 



c 


From this, by the elimination of p, is obtained 


P 


r 



2 


K-l 


'm 


Ec 1 /*-! 



The first factor is equal to the reservoir pressure p , as appears 
by setting q equal to zero. According to this 


K 



From the adiabatic equation it then follows that 

1 



where p Q is the reservoir density. 


Translation by Edward S. Shafer, 
National Advisory Committee 
for Aeronautics. 




